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Abstract 

Cosmological models with Lagrange multipliers are appealing because they could explain the behaviour 
of the dark sector in a unified way. In this work we analyse extensions to the "Dust of Dark Energy model" 
■ proposed in [l| by including spatial curvature and more general potentials of the scalar field. We perform 



in 



Of 



dynamical system analysis and we determine the evolution of the equation of state parameter as a function 
of the scale factor. We present observational constraints on this model by using Union2.1 dataset and H(z) 
data. 
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I. INTRODUCTION 



In 1998 measurements of the luminosity distance of supernovae type la (Snla) indicated the 
unexpected result that the Universe is undergoing accelerated expansion Q] , which would be driven 
by a negative-pressure matter component called dark energy. 

On the other hand, astrophysical observations provide compelling evidence for the existence 
of a non-baryonic, non-interacting and pressure-less component of the Universe, dubbed dark 
matter. This component clusters allowing structures to form. 

The existence of both dark components is supported by observations such as cosmic microwave 
background (CMB) 4\, baryon acoustic oscillations (BAO) Hubble constant measurements Q], 
Snla The available data indicate that 72.9% of the total matter content is dark energy, 

22.6% is dark matter and the remaining 4.5% corresponds to baryonic matter [4]. 

Given that we can not measure direct evidence of dark matter or dark energy, there exist a 
degeneracy in the dark sector 8(]. This degeneracy allow us to explore different dark candidates 
for the matter content of our universe. In this sense, it is very appealing to consider a single fluid 
describing the dark sectors in a unified way, which behaves as dark matter in early epochs and 
as a mixture of dark matter and dark energy nowadays. The archetypical unified model is the 
Chaplygin Gas ({J, which has been widely studied in several context [xo| . 

Recently, a new model for unifying the dark sector has been proposed This model, named 
Dust of Dark Energy (DDE), describes the dark sector by using two scalar fields where one of them 
is a Lagrange multiplier which imposes a constraint on the dynamics. In this sense, the dark sector 
is described by a single fluid which could represent dust or dark energy in different epochs of the 
evolution. The DDE model is appealing because it could be consistent with structure formation 



as suggested in Ref. 1 



. Cosmological models with Lagrange multipliers (LM) has been studied in 



different context [ll|- 15]. For example, in 2] the role of LM was analyse in the context of f(R) 
gravity, in [l| the Hamiltonian formalism was developed in modified gravity, in |l^| the authors 
investigate cyclic and singularity free scenarios in the context of modified gravity with LM and in 
[15I cosmological models with LM were studied with the focus on the cosmological constant value. 

On the other hand, the more recent cosmological data seem to favour a slightly closed geometry 
for our universe. A joint analysis with CMB, BAO and Snla indicate: ^ = -0.00571°;™ Q. 

In this work we extend the study of the model developed in [l| by considering: more general 
potentials for the scalar field and spatial curvature. We perform dynamical system analysis and 
put constraints on the parameters of the model by using the Union 2.1 sample of Snla and the 



2 



expansion rate data H(z). This paper is structured as follows: in section [XT] we describe the model, 
in section IIIII we use dynamical system analysis in order to study the asymptotic behaviour of 
the model. In section IIVI we show the numerical solution to the differential equations describing 
the model, in section [V] we perform Bayesian analysis with supernovae and H(z) data. Finally, in 
section [VI] we resume our results. 



II. THE MODEL 

The model is described by the action 

S = fd^x^-g [-A + K (V^) + a(X-^ 2 (V0) 

where K is a function of the scalar field ip and the kinetic term X = \g a ^V a ip^ pip- A is a Lagrange 
multiplier, p? is the potential of the scalar field ip, R is the Ricci scalar, g is the metric determinant, 
k 2 is a normalization constant and we consider c = 1. From this action we get the following set of 
field equations: 



G^ = k 2 T^ V,T; = and X = -p\i>) 



where is the Einstein tensor and = diag(p, — p, —p, —p) is a perfect fluid-type energy- 
momentum tensor. The total energy density p and the total pressure p are defined respectively 
given by [lj]: 

p = p 2 (K x + A) - K and p = K (1) 

In order to get the explicit form of the field equations we consider the Friedmann- Robertson- Walker 
metric in co-moving coordinates with a non zero curvature term: 

ds * = dt 2 _ a 2 {f) ( + r 2 dQ 2 + r 2 gin 2 

\ 1 — kr z 

where a(t) is the scale factor and the curvature parameter k = 0, +1, —1 represents flat, closed and 
open spatial sections, respectively. 

By imposing homogeneity and isotropy to the field equations we can consider X = ^ip 2 and, 
because of the constraint, p = ip. Dots denote derivatives respect to the cosmological time. 
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By choosing K = —X we recover the dynamics of the dark sector, where P = \ (2A — 1), 

2 

p = —\ and the derivative of the state parameter lj = - turns to be: 



Here primes denote a derivative respect to log a. We note that in order to have p > we need to 
fulfil the condition A > ^ or equivalently oj < 0, which will be assumed from now on. We recover 
a cosmological constant- type fluid for a constant potential p 2 = p^ an d A = 1, given that in this 
case p = 4p- and u) = — 1. 

By combining the conservation equation, p' + 3(p+p) = 0, the Friedmann equation, H 2 + -^ = | 
and the definition of p in Eq.([T]) we get: 

^3 



A' = (1 - 2A) 



—6u 
1 - kx 



(3) 



where we have defined e = — ^ and % = ~^r- denotes the derivative of respect to the scalar 
field ip and = - is the expansion rate. We have used k = 1. 

The dynamical set of equations describing the model is conveniently chosen to be given in terms 
of the functions co, \ an d e as: 
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CO 



3uj 



1 + w 



-6w 



X' = x(l + 3u/> 
e ' = _ e 2 (r _ 1} 



— 6cj 
1 - kx 



kx 



(4) 
(5) 
(6) 



where we defined T 



■4> 



We note that by assuming T as a function of e it is possible to get a 



closed set of equations. By providing T(e), the function p(ip) is determined by the solution of the 
following differential equation, T{e)p^ = p^p, where we have to use the definition of e in terms 
of the potential p. In this sense, to provide T(e) is equivalent to define the scalar field potential p. 
See TABLE [TT] for simple examples. 



III. DYNAMICAL SYSTEM ANALYSIS 



Constant e 



In order to study the set of Eqs.dU)-© we begin by considering the simplest case of a constant 
e. There are two possibilities to get a constant e: a potential p{ip) = Cte. which implies e = or 
p(ip) = Ae B ^ , for a non-zero constant e. In both cases the Eq.([6]) is trivially satisfied. 

Under these considerations the dynamical set of Eqs.@-([6]) is reduced to a two dimensional 
system described by: 



X 



3w 
(l + 3o;)x 




OJf y/l-kx 



(7) 
(8) 



where we have conveniently define ojj by e = ^^-^==, see [lj]. From this definition we note that 
for a constant potential ujf = — 1, otherwise cot will be a negative constant in order to have a real 
valued e. 

The critical points of the system and their main characteristics are given in TABLE U where the 
conditions for the existence of the critical points are shown. The most interesting critical points 
are 1 and 2 which could be consider as the past and future evolution of the universe, respectively. 
In this case the universe evolves from a state dominated by a fluid with u = (dust) to a state 
dominated by a fluid with to = tOf < —1/3 (dark-energy). This behaviour is corroborated by the 
numerical integration of Eqs.([7])-([8|), see FIGs. [T][3j 
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TABLE I: Critical points and stability behaviour for the system of Eqs. ([?])-©• We have considered \ > 
and ujf < 0, which follow from the definition of these variables. 

As an example, we show in FIGs. [T]l3]the phase space for four numerical solution to Eqs.([7])-([8]) 
and different values of the parameters uif and k. In these figures we have included the Direction 
Field of the system in order to have a picture of whatever a general solution looks like. In particular, 
in FIG. Q] it is shown the case where uif = — 1. In this case the curvature term is irrelevant, as 
we note from Eqs.([7])-([8]), and we reproduce the result in [lj]. Note that in the figure the physical 
part of the plot is delimited by u < and x ^ 0. In FIG. [2] the cases k = 1 for, ujj = —0.9 and 
ujf = —1.1 are shown. In FIG. [3]the cases k = —1 for us* = —0.9 and ujf = —1.1 are shown. 

It is interesting to note that for a nearly flat scalar potential where e<l and approximately 
constant, the system of Eqs.(jU)-([6]) can be reduced to a two dimensional system [la ]. The critical 
points consistent with this kind of potential are 1 and 2 of TABLE HI because in this case uif has 
to be close to —1. 

In section [V] we are going to contrast this particular model to cosmological observations via 
Bayesian methods, which will allow us to constraint the parameters of the model. 
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B. Variable e 



In order to study more general behaviour for the solution of Eqs.Q-©, where e is not a 
constant, we consider a family of potentials (see TABLE [TIj) which generate a simple structure for 
the (r — 1) term in Eq.© as (T — 1) oc e n where n is an integer and n > — 1. 

This family of potentials allows that the set of Eqs.([H)-([6]) becomes a three dimensional au- 
tonomous system with a critical point for e = 0, x = an d w = or w = -1. 



r - 


i 


= -n- 1 -> 




= /ioV'™ 




r - 
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= e- x -> 








r - 
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(( re + l)(V. + C 2 ))™/<™ + 1 ) 
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, ((n + l)(^-C 2 ))'*/(' l + 1 ' 
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TABLE II: The family of potentials which close the autonomous system of Eqs.(j4l)-([6j) 




FIG. 1: The Direction Field of Eqs. (H)-© and the numerical solution for 4 different sets of initial conditions. 
Case tjf = —1. 



As an example, let us show the case of a potential such that V — 1 
dynamical set of equations becomes: 



X. 



3w 1 + ui 

(l + 3w)x 
— e 



-6w 



1-kx 



6oj 



1-kx 



e 1 . In this case, the 
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FIG. 2: The Direction Field of Eqs. ©-© and the numerical solution for 4 different sets of initial conditions. 
Case k = 1, left panel w/ = —0.9 and right panel uif = —1.1 




FIG. 3: The Direction Field of Eqs. ©-© and the numerical solution for 4 different sets of initial conditions. 
Case k = —1, left panel ujf — —0.9 and right panel w/ = —1.1 



This system have two critical points given in Table ITO! The critical point 1* is an unstable focus 
and the critical point 2* is an attractor. Similar to the case discussed above, we can consider these 
critical points as the past and future evolution of the universe and both are independent of the 
curvature. In this case the universe evolves from a state dominated by a fluid with to = to a state 
dominated by a fluid with equation the state oj = — 1 (cosmological constant). This behaviour 
is corroborated by the numerical integration of Eqs.([9])- (|lip . Some of these solution are given in 
FIG. |U where a projection of the solution to the axes (w, x) an d (w, e) is shown, for several choices 
of the initial conditions. 
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TABLE III: Critical points and stability behaviour for the system of Eqs. ([§))- ([TTj) . 




1.0 -0.8 -0.6 -0.4 -0.2 0.0 



0.8 -0.6 -0.4 -0.2 



FIG. 4: Some numerical solutions for Eqs.([9|)- (|11[) where we have consider as an example k = —1 and different 
initial conditions. The left panel is a projection to the axis (w, \) an< i t ne right panel is a projection to the 
axis (lu, e). 

IV. NUMERICAL SOLUTION FOR A NEARLY FLAT SCALAR POTENTIAL 



In order to numerically integrate the set of Eqs.([7|)-([8|) we use more convenient functions defined 
and Qj- = — ai k H $ at any time in the evolution. In terms of these new variables Eqs.([7J)- 
are transformed to: 



as Q — 



3co 



i + u)- .[^-(i + uf)Vn 
V w / J 



n' = i)(i + 3w) 



(12) 
(13) 



By considering contributions to the spatial curvature of order £lk = ±0.005 (consistent with 
the data Q]) there are no significant modifications in the evolution of u{a) as it is shown in FIG. 
[SJ Here the subscript denotes the value of a function today. 

In FIG. [5] we see that the fluid behaved like dust in the past (a < ao), with a state parameter 
close to zero. This behaviour is independent of the allowed value of ujf and consistent with small 
contributions of spatial curvature, as we noted in the dynamical system analysis. In the future 
(a > ao) the state parameter reaches a constant value corresponding to u) = ujf, which is also 
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-0.645 




-2 -1 1 2 -0.210 - 0.208 - 0.206 - 0.204 - 0.202 - 0.200 

log[a/ao] log[a/a ] 



FIG. 5: The left panel corresponds to the evolution of state parameter for a flat spatial section. We have 
used w(0) = ujq = —0.74 which is consistent with the cosmological data that we present in section [V] The 
continuous line corresponds to ujj = —1.1, the dotted and dashed lines are for ujj — —1 and w/ = —0.9 
respectively. The right panel shows a closed region for uif = —0.9 and £lk = +0.005, 0, —0.005 from top to 
bottom. 

independent of the curvature. When ujf = — 1 the fluid asymptotically becomes a cosmological 
constant, whereas for cot < — 1 the fluid asymptotically becomes a phantom fluid. 

The curves in the left panel of FIG. [6] are the result of numerical integration of Eqs. (|12|) - ()13p 
with cj(0) = uq = —0.74 for different values of f2/% . The area between the solid lines expands a 
continuos range of values of the curvature parameter, inside the current observational constraints 
Q]. As we noted in the dynamical system analysis, the effect of curvature is not significant in the 
initial or final state of the universe in the context of this model. 

The right panel of FIG.[6]shows the numerical integration of the energy density p. For ojf = —1.1 
and ujq = —0.74, the final state of our universe will be dominated by a cosmological constant fluid. 
For — | < uif < —1 the energy density goes to zero at the end and for ojf < —1 we have a universe 
which will be dominated by a phantom fluid at the end, where the energy density p increases in 
the future without bound. 

V. OBSERVATIONAL DATA ANALYSIS 

In this section we examine the observational constraints on the model defined by Eqs.([7|)-([8]), 
with and without spatial curvature. We use Snla observations and H(z) data. 
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FIG. 6: The left panel shows the evolution of the density function f2 for three different values of J7fc . The 
continuous lines correspond to ojq = —0.74 and luj = —1.1, for Qj. = ±0.005. The dotdashed line is for 
ilk = 0- The right panel shows the evolution of the energy density with a;/ = —1.3, —1, —0.8 for solid, dotted 
and dashed lines, respectively. We normalized to have p = 1 today and we have considered ujq = —0.74. 

We perform bayesian statistical analysis using Snla data from the Supernova Cosmology Project 
Union2.1 sample [Tj], with 580 supernovae over the range 0.015 < z < 1.414. 
We fit the theoretical distance modulus jlth{ z ) defined by: 

'H d L (z)~ 



hh{z) = 5 log 



10 



+ Mo 



to the corresponding observed distance modulus fi bs,i- 

Here j2 = 42.38 - 51og 10 /i, H = W0h 
[kms~ 1 Mpc _1 ] is the Hubble constant and the luminosity distance is defined as cIl{z) = (1 + z)r(z) 



with 



r(z) 



-Sk 



n 



ka 



dz 



for Sk(x) = sin(x), x, sinh(x) for k > 0, k = 0, k < 0, respectively. 

The constraints from the Snla data can be obtained by minimizing the following x 2 function: 

580 , .. ~ , ~ , x 2 

M«.p) = E( ^5 J 

where p represents the model parameters and (Tfj,(zi) is the distance-modulus uncertainty for the 
corresponding redshift Z{. 

It is not difficult to realize that /jq is a nuisance parameter and we can easily marginalize over 



it Q. 



Thus instead of minimizing we minimize the function x~ which is independent of the 
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fio parameter. 
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FIG. 7: Contour plots (la and 2a) in the u>o — ujf and Clk — ujf plane for the joint constraint Snla+H(z). 
The best fit parameters are indicated with dashed lines. The gray region is excluded given the condition 
U)f < ujq. The left panel is for the flat case and in the right panel we have chosen wq = —0.736 as a prior. 
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Snla+H(z) 


570.974 


(prior) 


-0.736±0.061 


-1.038±0.193 


Snla+H(z) 


570.963 


-0.015±0.139 


-0.736 (prior) 


-1.040±0.062 



TABLE IV: The table shows the best fit parameters with the la uncertainty. We have performed Bayesian 
analysis with two free parameters and different prior in each case. In both cases we have imposed ujq > w/. 

We also perform statistical analysis using the Hubble expansion rate data {19I ]. In the same 
way as it was done with the jl^ parameter, we note that Hq is a nuisance paremeter and, instead 
of minimize the function: 

2 (tt \ ( Hobs,i - H t h(zi] H ;p) \ 2 
Xh(H , p) = > : — , where H th = H f(p) 

V <J H {Zi) ) 
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where now 
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&H(Zi) 2 
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; B(p) = £ 



c(p) = £ 



i=l 



i=l 



t=l 



A joint analysis using Snla+H(z) lead us to the best fit parameter showed in TABLE ll VI and FIGJ71 
VI. CONCLUSION 

We explore an alternative scheme for the problem of the dark sectors in Cosmology where dark 
energy and dark matter are described by the evolution of a single fluid. In particular, we analyse 
extensions to the DDE model proposed in jl| by including spatial curvature and more general 
potentials. 

We have found a family of potentials for which the model can be described by a three-dimensional 
autonomous system. We study the corresponding critical points and their characteristics. In 
general, there are two critical points which can be interpreted as the initial and final state of our 
universe. Namely, in the initial state the fluid behaved like dust whereas in the final state the 
fluid have a constant state parameter with ujf < — |. These results are independent of the spatial 
curvature. 

In order to constraint the parameters of the model by using Bayesian analysis, we find numerical 
solutions to the set of Eqs. (jl2fl - (|13|) . We found that the curvature has a negligible incidence in the 
evolution of the state parameter u, as it is shown in FIG. [SJ 

The Bayesian analysis shows that this model is consistent with the available data from Snla 
and H(z). For null spatial curvature the best fit values for the parameters are loq = — 0.736 ± 0.061 
and cjf = —1.038 ± 0.193, which suggest a final state dominated by a phantom-type fluid and a 
crossing of the so called phantom barrier in the future evolution. We note that the value of wq is 
consistent with the best fit results for ACDM when a single effective fluid is considered [4]. 

When spatial curvature is taken into account (choosing as prior ojq = —0.736), the best fit 
values are £lk = —0.015 ± 0.139 and ojf = —1.040 ± 0.062, slightly favouring a closed model. In 
both cases (with and without curvature), the best fit value for cjf is consistent with a nearly flat 
scalar potential. 

In the near future we expect to study cosmological perturbations on this model in order to 
explore possible deviations of the standard picture and include more data as WMAP and BAO. 
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